Submitted to A J 

Preprint typeset using 1^1^^ style eniulatoapj v. 5/2/11 



(N 

o 
a^ 

(N 



O 

u 

o 



"GALAXY," DEFINED 

B. Willman\ J.Strader^ 

Submitted to AJ 

ABSTRACT 

A growing number of low luminosity and low surface brightness astronomical objects challenge tra- 
ditional notions of both galaxies and star clusters. To address this challenge, we propose a definition 
of galaxy that does not depend on a cold dark matter model of the universe: A galaxy is a gravitation- 
ally bound collection of stars whose properties cannot be explained by a combination of baryons and 
Newton's laws of gravity. After exploring several possible observational diagnostics of this definition, 
we critically examine the classification of ultra-faint dwarfs, globular clusters, ultra-compact dwarfs, 
and tidal dwarfs. While kinematic studies provide an effective diagnostic of the definition in many 
regimes, they can be less useful for compact or very faint systems. To explore the utility of using the 
[Fe/H] spread as a complementary diagnostic, we use published spectroscopic [Fe/H] measurements 
of 16 Milky Way dwarfs and 24 globular clusters to uniformly calculate their [Fe/H] spreads and 
associated uncertainties. Our principal results are: (i) no known, old star cluster less luminous then 
My = —10 has a significant (J^O.l dex) spread in its iron abundance; (ii) known ultra-faint dwarf 
galaxies can be unambiguously classified with a combination of kinematic and [Fe/H] observations; 
(iii) the observed [Fe/H] spreads in massive (> lO^M©) globular clusters do not necessarily imply that 
they are the stripped nuclei of dwarfs, nor a need for dark matter; and (iv) if ultra-compact dwarf 
galaxies reside in dark matter halos akin to those of ultra-faint dwarfs of the same half-light radii, then 
they will show no clear dynamical signature of dark matter. We suggest several measurements that 
may assist the future classification of massive globular clusters, ultra-compact dwarfs, and ultra-faint 
galaxies. Our galaxy definition is designed to be independent of the details of current observations 
and models, while our proposed diagnostics can be refined or replaced as our understanding of the 
universe evolves. 
Subject headings: galaxies: star clusters — galaxies: dwarf — galaxies: kinematics and dynamics 



(N 
> 

oo 

o 

(N 

rn 
O 
(N 



% 



1. INTRODUCTION 

It h as been nearly a century since the term "galaxy" 
(|de Si tter 1917; Crommelin 1918; Shapley 1919) was first 
applied to the spiral nebulae that were later fou i id to be 
true extra-gala ctic stellar syste ms (|Slipheiill917l : iHubbld 
119261 see also iGrahamI 120111 and references therein). 
The term "cluster" has been used to refer to Milky 
Way (MW) open and globular cluster star syste ms since 
their initial discoveries more than 200 years ago (JMessier 
Il781h . In the early 20th century, the primary differences 
between objects classified as galaxies and star clusters 
were (i) the milkj0 appearances of galaxies versus the 
grainy appearances of star clusters, and (ii) the "island 
universe" environments of galaxies versus the association 
of star clusters with the MW system. In the intervening 
years, galaxies and star clusters have largely been clas- 
sified based on their physical sizes: galaxies have typical 
sizes of hundreds of pc to tens of kpc whereas star clus- 
ters have typical sizes of a few pc, with a scatter to tens 
of pc. 

The lion's share of known star clusters and galaxies 
can be classified by this simple "I know it when I see 
it" size-based distinction. However, there are a growing 
number of astronomical objects that are not so easily 
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•f The word "galaxy" derives directly from the Greek word for 
"milky" . 



classified — those at extreme low luminosities and surface 
brightnesses, and those f illing gaps previously o bserved 
in size-luminosity space (jMisgeld fc Hilkei)l2011fl . These 
objects currently hold the most promise to shed new light 
on galaxy formation at the bottom of the hierarchy and 
the distribution of dark matter to the smallest possible 
size and mass scales. For example, ultra-compact dwarfs 
(UCDs) have luminosities (—13 < My < —9) similar 
to those of both dwarf spheroidal galaxies and lumi- 
nous globular clusters (GCs), but sizes (10 pc < rhaif 
< 100 pc) intermediate to both populations and dynam- 
ical mass-to-light r atios of ^ 2 — gl, larger than those 
typic al of GCs (e.g..lHilker et al.lll999l:lDrinkwater et al.l 
120031: iHasegan et al.l 120051: lEvstimeeva et all 12007 



Mieske et al.li2008HChilingarian et aLll2011l:lBrodie et alj 
1201 1[ ). Some uhra- faint MW satehites (e.g.. Segue 1, 
Segue 2, Bootes II, Willman 1) have also challenged 
our notion of galaxies. These objects have luminosities 
(— 3 < Mv < — 1) lower than those of nearly any known 
old star cluster or dwarf galaxy, physical sizes (20 pc 
< rhaif < 40 pc) between those of most star clusters 
and dwarf gala xies, and dynamica l mass-to-light ratios 
as high as 3000 (jSimoiLet al."2011). While tidal dwarfed 
may provide a piece to this puzzle, they hayerooved dif- 
ficult to study and classify themselves (see lDudl2012l for 
a review). 

§ All mass-to-light ratios are in V unless otherwise stated. 

^ We use "tidal dwarf" rather than "tidal dwarf galaxy" 
throughout this paper, to avoid presupposing a galaxy definition 
for this class of objects. 
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The origin and properties of systems such as UCDs, 
extreme MW dwarf sateUites, and tidal dwarfs are fun- 
damental to many open questions in galaxy forma- 
tion and cosmology. They might hold unique clues 
to relationships between different classes of hot stel- 
lar systems (e.g., giant elliptical galaxies, dwarf ellip- 
tical galaxies, dwarf spher oidal galaxies, UCD s , nu- 
clear star clust e rs, GCs: iDabringhausen et al.l [20081: 
Wolf et al.l[20T0l: IMisgeld fc Hilkei<l20lH IZaritskv et al.l 
201 It ). They might be our best luminous tracers of 
sub-galactic dark matter. Having a well-defined clas- 
sification scheme will be essential to these studies, 
since imminent and u pcoming wide-field surveys, includ- 
ing Pan-ST ARRS 1 (iKaiser et al.l [200^ . the Southern 
Sky Survey fKeller et al."2 007D, the Dark E nergy Sur- 
vey ([The Dark Encrgv Surv ey Collaboration! 2005) and 
LSST (jlvezic et al.l I2008<) . are expected to reveal large 
numbers of previously unseen low surface-brightness sys- 
tems. 

As the rate of discoveries and the diversity of the 
known cosmic zoo increases, the question "What is 
a galaxy?" is bei ng discussed at conferences and in 
the literature (e.g.. iGilmore et al. 2007: Kr oupa 2008; 
Ivan den Be rgh 20081: iPb'rbes fc Kroupall20TTI ). The most 
common distinction currently made between galaxies and 
star clusters is the presence of dark matter — galaxies re- 
side at the cent ers of dark raatter halos and star clus- 
ters do not fe.g..lSimon et al.ll201ll: iTollerud et~aLll2011l: 



iWillman et al.l 12011 '). One strength of this definition is 
that it facilitates studies of dwarf galaxies in a cosmolog- 
ical context: it allows a straightforward connection be- 
tween the set of objects classified as galaxies and a com- 
parison with the predictions of dark matter plus galaxy 
formation models. One weakness of this physically mo- 
tivated definition is the fact that the cold dark matter 
model is a theory rather than a physical law. 

Other recently proposed definitions for a galaxy have 
included rhaif > 100 pc, a relaxation time longer 
than a Hubble time, o r complex stellar populations 
(jForbes fc Kroupal [20lTI ) . Ahhough each of these def- 
initions have their own strengths (namely that they 
are straightforward to diagnose), they each also have 
shortcoming!**!. For example, size-based classifications 
are becoming increasingly arbitrary as size-luminosity 
space is becoming cont inuously populated with objects 
(IMisgeld fc Hilkeiil2011[ ). The concept of "complex" stel- 
lar populations has also become ill-defined now that 
light element abundance spreads have been identified 
in a large frac tion of M W G C s (e.g., iGratton et al.! 
200l ICohen fc Mclcndei [2005 ID'Antona et al.! 120051: 



Carretta et al.! l2009b) . which were once thought to be 



pristine examples of simple stellar populations. 

In this paper, we tweak past definitions of the term 
"galaxy" with the aim of appealing to a sufficiently broad 
cross-section of astronomers that consensus might be 
reached. In §2, we motivate the importance of hav- 
ing a clear definition of galaxy within the astronomi- 
cal community, and then present a physically motivated 
definition. In §3, we consider kinematics and comple- 
mentary indirect diagnostics such as [Fe/H] spread and 
population-based diagnostics that might be used to test 

** Se e also the related, informal discussion linked at the URL 
|http : //arxiv . org/abs/1101 . 3309 . 



whether an object is a galaxy. In §4 we examine the 
known properties of ultra-faint dwarfs, UCDs, GCs, and 
tidal dwarfs in the context of our proposed definition. 

2. GALAXY DEFINITION 

The classification of astronomical objects is more than 
a semantic pursuit. Words matter. The terminology we 
choose to describe our research affects how appealing and 
accessible it sounds to funding agencies and the public. 
Words also dramatically affect the research community's 
ability to draw global conclusions from diverse sets of 
astronomical objects. For example, massive star clus- 
ter wCen may be the remnan t of a stripped, nu cleated 
dwarf galaxy (|Lee et al.l [19991: iBekki fc Freeman! [2003) . 
However, it is cataloged as a GC and thus is not con- 
sidered in studies of the MW's dwarf galaxy population. 
Well-defined, well-chosen classification schemes therefore 
improve our understanding of the universe by facilitat- 
ing meaningful comparisons between models and obser- 
vations. Conversely, ill-defined, ill-chosen classification 
schemes can muddy our understanding of astrophysical 
phenomena. 

Although no single definition of galaxy will be optimal 
for all purposes, we propose a physically motivated defi- 
nition that will facilitate studies of galaxies both in and 
out of a cosmological context: 

A galaxy is a gravitationally bound collection of stars 

whose properties cannot be explained by a combination 

of baryons and Newton's laws of gravity. 

In a dark matter context (whether cold, warm, self- 
interacting, or other) , this definition loosely translates to 
measuring whether an object contains dark matter. Al- 
ternatively, the definition can be interpreted to delineate 
those objects for whic h non-standard theories of gravity 
could be relevant fe.g.. lMilgromll983tlSotiriou fc Faraonil 
120101) . For such theories, our definition would require dis- 
tinct observable consequences of non-standard gravity to 
be imprinted on an astrophysical system for it to be clas- 
sified as a galaxy. In the interest of simplicity, we focus 
Sections 3 and 4 of this paper on diagnosing our pro- 
posed galaxy definition in a dark matter-based context. 
However, we encourage others to more explicitly explore 
these diagnostics in alternative contexts. 

We refer to Newton's laws on a macroscopic scale; this 
part of the definition should be considered to include ob- 
jects that require general relativity to be understood. We 
recognize that galaxies can have a significant amount of 
baryonic mass in non-stellar forms such as gas and dust. 
It is uncertain whether the local universe harbors any 
"dark galaxies" that contain gas but are entirely fre e 
of stars (Min chin etUI 120051: IDuc fc Bournaudi [2008). 
Whether these objects should be classified as galaxies 
is an open question that we do not confront here but 
that deserves further debate. 

A purely descriptive astronomical classification (such 
as relaxation time) may be relatively straightforward for 
either observers or theorists to implement. A weakness 
of the definition proposed here thus lies in finding ade- 
quate diagnostics to measure whether the properties of 
the lowest luminosity and most compact objects are ex- 
plicable with baryons and Newton's Laws (see §3 for a 
discussion of possible diagnostics) . 
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Because galaxy formation itself is hierarchical in a cold 
dark matter universe, there is no trivial distinction be- 
tween a single galaxy with satellite galaxies (such as the 
MW) and a galaxy group or cluster. As a diagnostic, 
IBusha et alj (J2012I ) propose that a bound collection of 
galaxies is a "galaxy" rather than a galaxy cluster if at 
least 50% of the stellar light is associated with one cen- 
tral object. This diagnostic is ultimately driven by the 
decreasing efficiency of galaxy formation in more massive 
dark matter h alos. 

The Busha et al.l (|2012f ) diagnostic helps guide intu- 
ition when exercising reasonable common sense in apply- 
ing our definition to astrophysical systems. For exam- 
ple, the intracluster star population of a galaxy cluster is 
composed of a gravitationally bound collection of stars 
whose dynamics cannot be explained by orbits within a 
Newtonian potential well dominated by cluster baryons. 
Such a system should not be classified as a galaxy be- 
cause it is physically associated with a galaxy cluster. 
Similarly, the Milky Way's stellar halo is a merely a com- 
ponent of the Milky Way — not its own galaxy. 

3. GALAXY DIAGNOSTICS 

In this and in the following section, we focus on objects 
for which a galaxy classification (or lack thereof) tends 
to be ambiguous. 

3.1. Stellar kinematics 

The most direct way to determine whether an object 
contains dark matter, or whether its properties are other- 
wise inconsistent with Newtonian gravity, is to conduct 
a kinematic study. The present day mass of a system 
is typically derived from its kinematics using formalism 
based on Newton's laws of gravity and the assumption of 
dynamical equilibrium. This dynamical mass can then be 
compared with the total mass present in the form of stars, 
stellar remnants, and gas. If dynamical mass exceeds the 
baryonic mass, then dark matter must be present or one 
of the dynamical assumptions — such as Newtonian grav- 
ity or virial equilibrium — must be fiawed. 

There are many regimes in which dynamical studies 
can be translate d with few assumptio ns into Newto - 
nian masse s (e.g. |WalkeLeL,al. 2009a; W olf et al.llMoh . 
IWolf et all ()2010h showed that the half-light mass of 
a dispersion supported system could be robustly cal- 
culated with only mild assumptions about the orbital 
anisotropy of its constituent stars. They derive M^aif = 
4G~^< afgg >rhaif- Here Mhaif is the total mass within 
the 3D deprojected half-light radius, < af^^ > is the 
luminosity weighted square of the line of sight velocity 
dispersion, and r^aif is the 2D projected half-light ra- 
dius. Such calculations have yielded (M/L )haif as high 
as ^3 000 for a MW satellite galaxy (Segue l. lSimon et al.l 

MM- 

It is not always possible to diagnose a galaxy defini- 
tion based on dynamical {M/L)haif alone. Many authors 
have looked at the relationship between M/L and other 
system properties (such as lu minosity, see e.g . Figures 3 
and 5 in lToUerud et al.|[20TTI and Figure 4 in 'Wo lf et~all 
l2010i ) While typical star clusters stand out as having low 
{M/L)haif (-- 1-5) for their luminosities (L ~ 10^-^Lq), 
dispersion supported galaxies {L ~ 1O*~^'^L0) have sim- 
ilar {M/ L)haif as star clusters. In such cases, a com- 



bination of {M / L)haif and other population arguments 
may be used to diagnose a galaxy classification (see also 
§3.3). Alternatively, dynamical modeling including trac- 
ers at larger distances can reveal M/L outside of r^aif- 
If the existence of dark matter is the correct inter- 
pretation of galaxy dynamics, then dynamical classifica- 
tion of galaxies may be robust to the effect of tidal mass 
loss. Simulations show that galaxies tidally stripped of 
mass should maint ain a high dynamical ma ss-to-light ra- 
tio. For example, i Penarrubia et all (|2008| ) showed that 
the mass-to-light ratios of tidally evolving dwarf galax- 
ies increase over time, assuming they reside in cuspy 
dark matter halos. Even if the dark matter halos host- 
ing dwarf galaxies are cored, their central dark matter 
densi ty slopes remain constant during tidal evolution 
(Penarrubia et al. 2010). 

3.1.1. Special Considerations 

Generally, dynamical M/L > 5 may be taken to diag- 
nose a galaxy classification, because such M/L are larger 
than expected from typical stellar populations. How- 
ever, a number of challenges face attempts to determine 
whether an observed dynamical AI/L of a system is con- 
sistent with expectations from baryons alone - especially 
for systems with M/L < 10, low intrinsic velocity disper- 
sions, or low surface brightness. Some of these challenges 
are fairly obvious, as the dynamical AI/L expected for a 
purely baryonic population varies significantly with: age, 
metallicity, initial mass function, dynamical state, and 
gas content. In this section, we highlight several specific 
examples which are less commonly discussed in the lit- 
erature. See also §4.2.1 for a more nuanced discussion of 
dynamical M/L in the context of UGDs. 

Several effects can cause an overestimate in dynam- 
ical mass, and thus an overestimate of M/L. For ex- 
ample, the orbital motions of binary stars can inflate a 
system's observed velocity dispersion. A recent, multi- 
epoch velocity study of Segue 1 suggests that binaries 
should not pose a major problem for the dynamical clas- 
sification of systems with in trinsic velocity dispersions of 
at lea st a few km s~^( Martinez et ani201lHSimon et al.l 
120 lit ). However, binaries d o materially impact lower 
velocity dispersion systems ([Bradford et al.|[2011[ ). and 
models based on more extreme assumptions than pre- 
viously considered identify regions of parameter space 
where binaries c ould impact Segue 1-likc velocity dis- 
persion systems (jMcConnachie Sz Co te 2010). Tidally 
unbound and MW foreground stars can also contami- 
nate spectroscopic samples of a MW companion's stars 
and inflate its observed velocity dispersion. The effect of 
such contaminants can be mitigated by a combination of 
careful simulations of the MW for eground and its color - 
magnitude-velocity distribution (jWillman et al.l 120 lit ) , 
the use of spectroscopic abundance indicators, statis- 
tical approaches to identifying object members (e.g., 
iWalker et al.l l2009b[) . and approaches to eliminating 
tidally strip ped stars that have been iii formed by N-body 
simulations (jKlimentowski et al.ll2007l ). 

Other effects may alternatively cause an underesti- 
mate of stellar mass, and thus an overestir nate of the 
presen ce of non-stellar mass. For example, iHernanded 
()2012D shows that the lowest luminosity systems (L ~ 
500 L0) can have {M/L)steiiar between 1 and 10 sim- 
ply from the stochastic effects of sampling an IMF with 
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a small number of stars. A tidally stripped, dynami- 
cally relaxed (and therefore mass-segregated) GC can 
also have a super-stellar M/L once the majority of its 
mass has been lost. Models of star cluster evolution 
that include the effects of mass segregation and the 
Galaxy's tidal field have shown that high fractions of 
stellar remnants accu mulate in the center as a clus- 
ter is stripped of mass (iVesperini fc Heggiil 19971 : iGierszl 
120011: iBaumgardt fc Makinoll2003[ ). Although possible, 
it should be rare to observe a system so tidally stripped 
that its global M/L is significantly inflated by this mech- 
anism. For example, although Palomar 5 is estimated 
to be ~100 Myr from complete destruction (less than 
1% o f its total lifetime), it is observed to have M /Ldyn 

< 1 (jOdenkirchen et all [20021: iDehnen et al|[200l . Ob- 
servational limitations may also generate ambiguity in 
the dynamical classification of the lowest luminosity {L 

< 1000 Lq) and low velocity dispersion (<3 km s~^) 
systems. For example. Segue 3 (L = 90^4° Lq, d 
'^ 17 kpc) contains only a few dozen member stars 
brighter than r = 22. 32 of Segue 3's stars were 
observed with Keck/DEIMOS to obtain velocity mea- 
surements with u ncertainties per exp osure per star of 
-3-10 km s-i ()Fadelvet al.|[20Tlf ). With a aios of 
0.3 km s~^ expected based on stars and Newtownian 
gravity alone, its measured velocity dispersion of 1.2 ± 
2.6 km s~^ is dynamically consistent with either a galaxy 
or a star cluster interpretation. Even with techniques 
which retrieve stellar velocities from me dium-resolution 

fsctra with uncertainties <1 km s~^ ()Koposov et aD 
11), star-poor systems need to reside within ^^20 kpc 
for there to be a sufficient number of stars bright enough 
to spectroscopically observe with high S/N with a 10m- 
class telescope. 

3.2. [Fe/H] Spread 

Another way to directly constrain the potential well 
in which a system formed is the presence of an [Fe/H] 
spread. The use of [Fe/H] as a diagnostic for our 
proposed galaxy definition is supported by a combi- 
nation of models of supernova winds in low-mass sys- 
tems and the observed abundances of stars in well- 
studied dwarfs and GCs. Iron is produced by super- 
novae (both Type II and la), so a dispersion in [Fe/H] 
implies that the system was able to retain supernova 
cjccta to form multiple generations of stars. The en- 
ergetic winds of supernovae can only be retained in 
a gravitational well of sufficient depth. Estimates for 
the GC r uass needed to retain SN ejecta are >fewxlO ^ 
Mq f^g.- lDopita fc Smith|[l98llBaumgardt et al.ll2008D . 
Observed [Fe/H] spreads of over 1 dex combined with 
inferred stellar masses of ~1000 Mq or less have thus 
contributed to a galaxy classi fication for both Se gue 1 
and Willman 1 (iMartin et a l. 2007^ jNorris et al.ll2010l : 
ISimon et al.ll201lHWillman"et al...20m 

3.2.1. Calculating crpe/H] 

To empirically investigate the difference in [Fe/H] 
spread, (T[fc/h] j between well-studied dwarf galaxies and 
GCs, we estimate the spread and associated uncertainty 
for each of 16 dwarfs and 24 GCs with publicly avail- 
able, spectroscopic [Fe/H] measurements. We only used 
[Fe/H] measurements based on actual iron lines, rather 



than studies that infer iron abundance from the calcium 
triplet or photometry. We used Bayesian Markov Chain 
Monte Carlo techniques to fit a normal distribution to the 
stellar [Fe/H] values for each object, taking into account 
the reported measurement uncertainties and assuming 
flat priorsO We summarize the standard deviation of 
each sample, (T[fo/h]i as the median of its posterior dis- 
tribution, together with a 68.2% credible interval (anal- 
ogous to the usual la confidence interval). Calculated 
values of [Fe/H], cr[Fc/H]i associated uncertainties, and 
references are summarized in Table 1. The uncertainties 
on the variances are an increasing function of decreasing 
sample size, because small samples poorly sample the 
underlying [Fe/H] distribution. 

A few notes on un usual cases: For S egue 1, we in- 
cluded the star from ISimon et al.l ()2011[ ) with only an 
upper limit to its [Fe/H] as a censored datum in our 
analysis. We used the largest set of [F e/H] values for 
wCentauri ([Johnson fc Pilachowskil[2010( ). However, this 
sample is biased against the most metal-rich subpopu- 
lation because it is magnitude-limited in V. We thus 
consider our estimate of i ts [Fe/ H] spread to be a lower 
hmit. The IMarino et all ([20111 ) data for M22 does not 
contain uncertainties, and so our reported cr[Fe/Hl is an 
upper limit. Our analysis does not include the Terzan 
5 GC despite claims of an [Fe/H] sprea d in this object 
([Ferraro et al.ll2009t lOriglia et aLll2011[ ). owing to its - 
solar abundance (and thus different origin than the old 
metal-poor stellar populations we are primarily consider- 
ing) and the possibility that the sample may be partially 
contaminated by bu lge stars. We also d id not include 
NGC 5824, in which iSaviane et al.l ([2012D have reported 
""[Fo/H] ^ 0-11 ~ 0-14 dex, because this measurement is 
based on the Calcium triplet (thus revealing a Ca spread, 
not necessarily an Fe spread). The GC NGC 2419 is 
kno wn to display a — 0.2 dex spread in Ca, but none in 
Fe ([Cohen etal. l [2?)T(il) . 

Although reasonable indicators of the dispersion in 
[Fc/H], the values in Table 1 should be considered with 
caution before comparing in detail with models. The ac- 
curacy of our estimates of the variance of [Fe/H] (and 
its uncertainty) rely on (i) the appropriateness of the 
underlying Gaussian model, (ii) clean membership sam- 
ples, and (iii) accurate uncertainties for individual stars. 
For the fai nter dwarfs in this set the first condition rarely 
holds (^e.g. JKirbv et al.|[201lD . so our estimated variances 
should be taken as indicators of the spread in metallic- 
ity rather than as exact values. The faintest dwarfs may 
also have a small number of interloper stars in their mem- 
bership samples (see also §3.2.4). The third condition — 
estimating accurate uncertainties — is most relevant for 
GCs, because their measured (Trpc/H] sue comparable to 
(or less than) the measurement uncertainties for single 
stars. For this paper, we have included the random un- 
certainty in the Fe I abundance a s the st andard error 
of the mean, while iCarretta et al.l ([2009a[ ) included no 
measurement uncertainties in their calculation of [Fe/H] 
spread. The practical effect is that the intrinsic [Fe/H] 
spreads we derive for GCs in this paper are slightly 
smaller than those in Carretta et al., (2009a), by typ- 
Carretta et al.l (|2009al) . we em- 



ically 0.01 dex. Like 



tt The relevant (simple) code is available on request. 
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phasize that their and our values are upper limits to 
be true [Fe/H] spreads because of our limited ability to 
mod el the full measureme nt uncertainties on each star 
(see iCarretta et al.l l2009al for detailed discussion of the 
relevant modeling issues). 

3.2.2. cr[Fo/H] m Mv > -10 Objects 

Figure [1] shows crpc/H] for MW dwarf galaxies (filled, 
black circles) and MW GCs (open red circles) as a func- 
tion of absolute magnitude. Uncertainty bars show the 
68.2% confidence intervals. We show objects with dy- 
namical classifications of galaxy or star cluster as differ- 
ent symbols in Figure [1] to highlight the regime in which 
"■[Fe/H] results in the same inference about a system's 
potential well as a dynamical study. This tests whether 
cr[Fo/H] niay be used as a galaxy diagnostic in cases where 
dynamical studies are inconclusive. 

This figure shows a striking difference between the 
[Fe/H] spreads observed for dwarf galaxies and GCs. The 
dwarfs all have spreads of 0.3-0.7 dex (even higher for 
Segue 1), whereas none of the GCs less luminous than 
My = — 10 have substantial [Fe/H] dispersions. After 
the upper limit of crpc/H] =0.1 dex estimated for M22, 
the next highest spread is 0.08 dex estimated for NGC 
6441. Although these values are small, they are formally 
greater than zero with > 99% probability (as calculated 
above) . These estimates may reflect the detection of mi- 
nor star-to-star variations in [Fc/H] in GCs less luminous 
than My = —10. However, in light of the caveats given 
above, they may yet be found to be consistent with no 
star-to-star variation in [Fe/H] . 

For objects less luminous than My = —10, the di- 
chotomy between (T[fc/h] of dwarf galaxies and GCs un- 
derscores that the dwarf galaxies formed within much 
deeper potential wells than the GCs. We conclude that 
a cr[Fo/H] > 0.2 dex in such systems would be sufficient to 
diagnose a galaxy classification because it would not be 
explicable with a combination of baryons and Newtonian 
gravity (without invoking substantial mass loss). While 
iron is not the only element that may provide relevant in- 
sight to the gravitational potential wells of objects in this 
luminosity regime, it is clear that iron spread provides a 
powerful diagnostic of the provenance of such objects. 

3.2.3. o-[Fe/H] m Mv < -10 Objects 

The interpretation of the iT[fc/h] spreads observed in 
the two GCs more luminous than My = —10, M54 
and wCen, is less straightforward. One interpretation 
of the spreads in M54 and wCen is that they are the nu- 
clear star cluster cores remaining from a stripped dwarf 
galaxy (M 54: S gr core,| Saraiedi ni fc LavdenlllQQSt wCen, 
iLee et aLlUggllBekki fc Frccmanl |2003[) . It remains to 
be seen whether the properties of the gravitationally 
bound remains of such a stripped galaxy would satisfy 
our definition of a galaxy, and be formally classified as 
such. Recent observations have discovered a significant 
amount of tidal debris that may be associated with wCen 
(jMaiewski et al.l[2012l ). Sgr is already a classified galaxy, 
so M54 would not be considered a separate entity. 

An alternative interpretation of the [Fe/H] spreads in 
these Mstar > 10^ Mq clusters is self-enrichment by 
SNe without the additional gravitational inffuence of 
dark matter or a non-Newtonian effect. This interpre- 
tation is complicated by the fact that M54 and wCen 



do not actually have the highest escape velocities of 
the GCs in our sample. Using the fitted relation be- 
tween central velocity dispersion and central esc ape ve- 
locity, Vesc,o/'^o = 3.7 + 0.9(c — 1.4), from.Gnedi n et al.l 
(j2002f ). we 'find that 8 of the 62 GCs with vel ocity disper- 
sions reported in the 2010 edition of iHarrisI ([1996!' ) have 
central escape velocities larger than M54's Vesc,o ~ 45 
km s~^(not including wCen). 5 of these (47 Tuc, NGC 
2808, NGC 6388, NGC 6441, and M15) are in our sample 
and do not display [Fe/H] spreads > 0.1 dex. NGC 6441 
and 6388 have escape velocities of 72 and 76 km s^^, 
respectively, larger than wCen's escape velocity of 61 
km s""'^. A caveat to this analysis is that these values are 
measured at the present day. At earlier times, these GCs 
were all more massive but have since undergone stellar 
evaporation and tidal mass loss; some may have also had 
different sizes. All of these factors could have affected 
their relative abilities to retain supernova ejecta. 

Observations of GCs in other galaxies provide tenta- 
tive support for self-enrichment in iron in My < —10 
GCs. For example, HST/ACS photometry of three of 
the most massive GCs in M31 are sugge stive of spreads 
in [Fe /H] on the red giant branch (|Fuentes-Carrera et al.l 
|2008() . The dynamical masses of these GCs range from 
2 — 6 X lO^Mp, co mparable to or larger than w Cen 
(IStrader et al.l[20ll . The M31 cluster Gl (3 x lO^M^;,) 
also may have a significant [Fe/H] spread (jMevlan et al.l 
I2OOI . 

Separately, a number of groups have identified evidence 
of self-enrichment in extragalactic GCs. Precise pho- 
tom etry of blue, metal-poor GCs in a variety of galax- 
jes (Harris et al. 2006; Mieske et al. 2006; Strade r et al 
2006; Spitler et al. 2006 ; Forbes et al. 2010 ; Mieske et al 



2010( ) shows a correlation between magnitude and color 
for metal-poor GCs. This mass-metallicity relationship 
is not observed in all galaxies studied, but a typical re- 
lation is Z ^ M"**, where Z is the mean metallicity of 
the GC and M is its mass. The onset of the correlation 
appears to be between ~ 2 x 10^ and 10^ Mq. The slope 
and onset mass of the correlation can be reasonably ex- 
plained by models in which the GCs self-enr ich in iron 
(|Bailin fc Hamsll2009HStrader fc Smithll2008n . 

If (nearly) all GCs with stellar masses above few 
xlO^ M0 display [Fe/H] spreads, then it is likely these 
spreads accrue from self-enrichment without the help of 
an additional gravitational field. More extensive spectro- 
scopic and photometric campaigns to quantify the [Fe/H] 
spreads of extragalactic GCs will be essential to develop 
a fuller picture of the connection between cr[Fe/H] and the 
formation channel(s) of objects with Afy < —10. 

3.2.4. A Relationship Between cr[p(,/Hi t^nd My For Dwarfs? 

Figure [1] displays another striking trend in addition 
to the dwarf/GC dichotomy: the apparent increase 
in g[Fe/H] w ith decreasing luminosity (see also §6.2 of 
iKirbvet al] 120 lit ). While the dispersion in [Fe/H] for 
most MW dwarf galaxies with My < —8 (the classi- 
cal dwarfs) is 0.3 — 0.4 dex, the dispersion for most of 
the lower luminosity dwarfs (the ultra-faint dwarfs) is 
0.5—0.6 dex. The most likely explanations for this appar- 
ent trend are: (i) a true physical difference in the crrpc/Hi 
of the least luminous systems, (ii) a systematic bias in 
the calculated ctjfo/h] as the model assumptions become 
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Fig. 1. — The dispersion in [Fe/H] measured for MW dwarf galaxies (black, filled) and globular clusters (red, open), calculated assuming 
an underlying Gaussian distribution. The systems shown with dwarf galaxy symbols in this figure have been independently classified as 
galaxies by dynamical studies. Willman 1 does not have a definitive dynamical classification, and so is shown as an open hexagon with a 
cross. By this figure, a galaxy classification can be indirectly inferred from Willman I's spread in [Fe/H]. The presence of a spread in [Fe/H] 
can diagnose a galaxy definition because it constrains the depth of the potential well in which a system formed, as supernova ejecta must be 
retained to form further generations of stars. Error bars show the estimated uncertainty on each dispersion given the [Fe/H ] measu rement 
uncertainties on the individual member stars. Values and references are summarized in Table 1. Figure 7 of lCarretta et aT] I l2010al ') shows 
a figure similar in spirit to this, but for a smaller set of objects and without measurement uncertainties. 

increasingly poor with decreasing luminosity, or (iii) a 
result of a low level of foreground contamination that 
disproportionately affects spectroscopic samples of the 
lowest surface brightness systems. The faintest dwarfs 
have tails at the metal-richer ends of their metallicity 
distribution functions that are not present in the classi- 
cal dwarfs. It is not yet clear whether those metal-richer 
tails are physical or a result of mild contamination in the 
samples. Exploring the relative likelihood of these three 
scenarios is beyond the scope of this paper, but will be 
imperative to pursue in the future. 



3.3. Indirect Diagnostics: Population Arguments 

Population arguments rely on the assumption of a sin- 
gle classification for all astrophysical objects known to 
populate a particular region of parameter space. Such 
arguments are handy because, for example, it would be 
both impractical and unnecessary to conduct a detailed 
analysis of each of the 200 million galaxies cataloged by 



the e ighth Sloan Digital Sky Survey data release ([Aiharal 
1201 If ) before classifying them as such. The most com- 
mon population-based classification is simply the size- 
based classification that is naturally made for galaxies 
with scale size > 1 kpc. All objects satisfying this "I 
know it when I see it" criterion that have been studied 
in sufficient detail have been kinematically shown to sat- 
isfy our proposed definition of galaxy (not including tidal 
dwarfs, see §4.4.) Some kinematic studies of galaxies 
have postulated that no unseen matter or modification 
of Newtonian gravity may be needed to explain their dy- 
namics fe.g.. iRomanowskv et al.|[2003D . However, such 
studies have alway s been shown to b e flawed on theoret- 
ical grounds (e.g.. iDekel et al.ll2005D or were refuted by 
subsequent observational studies. 

Attempts have been made to connect, or distinguish, 
galaxies and star clusters using scaling relations that 
combine their metallicities, effective sizes, internal ve- 
locities, luminosities or derivatives thereof. Such studies 
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have recently focused on variants of the Fundamental 
Plane such as the Fundam ental Manifold (jForbes et alj 
12011'; 'Zaritskv et"alll2011[ ) and the Fundamental Curve 
|Tollcrud ct al.l l2011| ). These scalings reveal similarities 
and differences in the ways baryons coalesce within dif- 
ferent types of systems. However, the scalings do not 
seem to shed light on the classification of objects as a star 
cluster or a dwarf gala xy, in a way niore meaningful than 
M/i. within r/,,„7f (e.^.jForbes et al.l2008HTollerud et alJ 



1201 It iZaritsk v et alJ 1201 ll ). One simple difference be 
tween galaxies and globular clusters as a population is 
the metallicity-luminosity relation observed for galaxies 
(but not star clusters, UCDs, or nearby tidal dwarfs) 



over a wide range of stellar masses (e.g. S killma n et al 
19891 iTremonti et"al|[20p: iWoo et al.|[2008t iKirbv et al 



201 It ) . Although the metallicity and luminosity of an in- 
dividual object would not be sufficient to classify it as a 
galaxy or star cluster, consideration of the metallicities 
and luminosities of a population of objects may aid in 
its classification (see also §4.2.2). It is also worthwhile to 
consider the placement of individual ambiguous objects 
with respect to observed scaling relations. Inconsistency 
with well-established relationships on a case-by-case ba- 
sis may be a sign that some of the cautions raised in 
§3.1.1 are affecting the kinematics, effective mass, or size 
measured for an object. 

Another approach to population-based classification is 
to include a broad set of properties such as spatial dis- 
tribution, metallicity, and orbits when looking for subtle 
trends within a diverse set of observables. The combi- 
nation of such a set of clues may help reveal whether 
some object or type of object with an ambiguous classi- 
fication has an origin (and thus, possibly, classification) 
mor e similar to t hat of star clusters or of dwarf galax- 
ies. iBrodie et al.l (j2011) recently conducted a thorough 
analysis of UCDs around M87 in the Virgo cluster. They 
combined size-luminosity, age-metallicity, spatial distri- 
bution, and orbital dynamics to infer the possible co- 
existence in size and luminosity of three sub-populations 
of UCDs: the stripped nuclei of dEs, remnants from more 
massive red galaxies (either their nuclei or merged clus- 
ters), and genuine star clusters. 

Although we do not aim to be exhaustive, throughout 
§4 we will mention some specific indirect diagnostics that 
may contribute to a galaxy classification. 

4. SOME EXAMPLES THROUGHOUT THE COSMIC ZOO 

In this section, we use the diagnostics in §3 to consider 
the classification of four populations of astrophysical ob- 
jects: extreme ultra-faint dwarf galaxies, UCDs, GCs, 
and tidal dwarfs. 

4.1. Ultra-faint dwarfs with ruaif < 50 pc 

We begin our discussion with extreme ultra-faint 
dwarfs, because their classification is starting to con- 
verge in the literature. The term "ultra-faint dwarfs" 
refers to the dwarf galaxies with absolute magnitudes 
fainter than My ^ —8. Currently, such objects are only 
known around the Milky Way and M31 because they are 
difficult to detect, although the Next Generation Virgo 
Cluster Survey should soon reveal them in Virgo. The 
most extreme of these objects (Segue 1, Segue 2, Bootes 
II, Willman 1) are observed to have My ^ —2.5 and 
rhaif ~ 30 pc. 



These extreme objects have total luminosities less than 
individual bright red giant branch stars. Their sizes are 
intermediate between typical GCs and low luminosity 
dwarf spheroidal galaxies. Despite their extreme and un- 
usual properties, direct and/or indirect diagnostics sup- 
port a galaxy definition for all four of these objects. With 
an {M/L)haif -- 3400 and cr[Fo/H] = 0.75J:|^;23 dex (Ta- 
ble 1), Segue 1 is a galaxy as diagnosed by both its kine- 
matic and its (T[fc/h]- Taken at face value, the dynam- 
ics of Willman I's stars require a high dynamical mass 
relative to its stellar mass. However, its irregular kine- 
matic distribution hinde rs drawing a robust cl assification 
from kinematics alone () Willman et al.l 120 lit) . Regard- 
less, the substantial spread in [Fe/H] among three mem- 
ber stars in Willman 1 (0.56^q 23 '^^^^ Table 1) demon- 
strates its galaxy classification. A dynar nical study based 
on sm all numbers of stars in Segue 2 (jBelokurov et al.l 
l2009f ) is consistent with a galaxy classification, although 
the uncertainties are still large. Finally, tidal arguments 
for Bootes II have suggested that it may need a sub- 
stantial dark matt er component for it to be self-bound 
(IWalshet al.l[2008l) . 

It is essential that all (candidate) extreme ultra-faint 
dwarfs close enough to study with lOm-class telescopes 
are spectroscopically investigated. Surveys like DES and 
LSST have the potential to uncover large numbers of 
objects like Segue 1 to distances beyond the reach of 
today's spectroscopic resources. If a sufficient number 
of nearby Segue 1-like objects are demonstrated to be 
galaxies, then systems discovered to share that region of 
size-luminosity space in the future might be classified 
as galaxies without extensive follow-up. Even now, it is 
not yet certain whether Segue 2 and Bootes II should 
be counted as galaxies or remnants thereof. Their clas- 
sifications will greatly impact the predicted number of 
luminous dwarfs orbiting the Milky Way and our (cur- 
rently minimal) knowledge of the bottom of the galaxy 
luminosity function. 

As a population, the MW's ultra-faint dwarfs follow 
lumi nosity- metallicity and luminosity- (M/L) relations 
(e.g.. lGeha et al.|[2009D . These scaling relations rule out 
pathological explanations for the ultra-faints as a popu- 
lation, such as clumps in tidal streams or stellar streams 
at orbital apocenter. When in doubt for any particular 
object, hypothesis testing against the spatial-kinematic 
predictions of a specific model can be use d to effectively 
vet a galaxy classification. For example, iZolotov et al.l 
(|2011| ) showed that the highly elliptical Hercules dwarf 
spheroidal is inconsistent with a cusp catastrophe hy- 
pothesis. 

4.2. UCDs 

Like "ultra-faint dwarf", the term UCD has no for- 
mal definition. It is generally used to refer to sys- 
tems with — 13 < Mv < —9 and 10 pc < rtaif < 
100 pc. This population of objects has proved partic- 
ularly challenging to classi fy. With up to 100 UCDs pos- 
sibly orbiting M87 alone (jBrodie et al.l 12011) . whether 
or not these should be counted as galaxies bears great 
importance for understanding the dwarf galaxy popu- 
lation of the Virgo Cluster in a cosmological context. 
Thus far, studies seem to be converging on the con- 
clusion that multiple formation channels may be re- 
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quired to explain the UCDs as a population, such as 
very massive s tar clusters or as the stripped nuclei of 
dwarf gala xies (iBrodie et al.ll2011l:lChiboucas et al.ll2011l : 
IChilingari an et alJl20lH iDa Rocha et all 1201 Ih . In this 
section, we do not review the work relying on popula- 
tion arguments or detailed s tudies of individual UCDs 
(e.g.. iMaraston et all 120041: iFellhauer fc Kroupal 120051: 
iNorris fc Kannappanll2011l ) to reach this conclusion. We 
instead discuss the efficacy of UCD kinematic studies in 
a cosmological context and consider possible future kine- 
matic and [Fe/H] UCD classification diagnostics. 

4.2.1. UCD Kinematics 

Dynamical studies of UCDs do not provide a clear 
diagnosis of a galaxy classification. UCDS, the most 
luminous UCD in the Fornax cluster (My — -13.55, 
fhaij = 87 pc), is the only UC D with spatially resolved 
kinematics (|Frank et al.l 120111 ). UCDS has less than a 
SS% mass contribution from dark matter within 200 pc, 
and M /L = S.6±0.S, if it is assumed that mass follows 
light. This M/L ma y be consistent with the M/L of 
its stellar popu lation fChilingarian ct al."2011', however 
see IMieske et a l. 2006 and Firth ct al. 2009 who esti- 
mate a lower stellar M/L). The spatially unresolved dy- 
namical studies of other UCDs yield dynamical M/L = 
2 — 5, plausibly (b ut not certainly) consistent with their 
stella r M/L (e.g., iHilker et all 119991: iDrinkwater et ahl 
20031: iHasegan et al.l 120051: lEvstigneeva et al.l 120071: 



Mieske et all 120081: IChilingarian et al.l 120111 ). The dy 



namical M/L of Virgo cluster UCDs seem to be system- 
atically higher than Fornax cluster UCDs (jMieske et al.l 
[2008,) ■ The inflated Virgo UCD M/Ls may be explained 
by unusual IMFs; a top-heav y IMF can yield a large frac- 
tion of dark stellar remnants ([Dabringhausen et al.ll200"9l 
[2012), while a bottom-heavy IMF is rich in low-mass M 
dwarfs with high individual M/L. 

It is not surprising that dynamical studies of UCDs 
do not easily yield a galaxy classification, even if (for 
example) they do presently reside in dark matter halos. 
To quantify this, we must begin with a reasonable hy- 
pothesis for the amount of dark matter expected within 
the half-light radii of UCDs if they reside in dark mat- 
ter halos. There are no simulations of sufficient spa- 
tial resolution to predict the expected amount of dark 
matter in the innermost ~S0 pc of a dark matter halo, 
and the highest resolution simulations do not include the 
effect of baryons, star formation, or feedback. More- 
over, there are known differences between the central 
mass densities observed for dwarf galaxies and the cen- 
tral dark matter densities predicte d for dwarf galaxies us- 
ing dark ma tter only simulations (jBovlan-Kolchin et al.l 
120111 I2012D that have, in son ie cases, been resolved 
with baryonic physics (e .g., iGoverna to ct al." 120101: 
iPontzen fc Governato 20121 ) . Similarly, Tollcrud et al.l 
(j2011| ) find that the observed central mass densities of 
UCDs are not consistent with residing in the Navarro- 
Frenk- White profile dark matter halos predicted by dark 
matter simulations. 

We therefore rely on an empirical hypothesis for the 
possible dark matter content of UCDs: they contain 
the same amount of dark matter within their half- 
light radii as known dwarf galaxies with the same half- 
light radii. We consider two MW dwarfs: Segue 1 
([Martin et all [20081: iSimon et al.l [20lll aios = ^-itW 



km s"\ My = -I.SIq 8, rhaif^ 29 ib 6 pc) and C oma 
Berenices ([Simon fc Gehal 120071 : [Mufioz et all 120101 My 
= _3_6±0.6, rhaif= 74 ± 4 pc). Using the I Wolf et"al] 
([20101 ) formula, we calculate Mseguei,haif = 3-7^2 J x 10^ 
Mq and McomBerMaif = 1.5±0.5x 10<^ M q . To obtain 
the half-light dark matter masses of these objects, we 
simply subtract out their approximate stellar masses as- 
suming a stellar AI/L of 2. Because Segue 1 and Coma 
Berenices are highly dark matter dominated, the derived 
dark matter masses depend little on our assumed value 

of {M/L)star- 

We use the half-light dark matter masses of Segue 1 
and Coma Berenices to predict the possible dynamical 
{M/L)haif of UCD- luminosity systems with half-light 
radii of 30 pc and 75 pc. Figure[2]shows the resulting pre- 
dictions, as a function of absolute magnitude and assum- 
ing a stellar M/L = 2 for the UCDs. UCDs are typically 
observed to have —13 < My < —9 and 10 pc < rhaif 
< 100 PC (see e.g.lMadr id etal.l[20Tol: IBrodie et al[ l201ll: 
iMisgeld et al.l20lH ). We predict that UCDs in dark halos 
would have dynamical M/L within their half-light radii 
of 2 — 3, consistent with observations. Given the large un- 
certainties in deriving stellar M/L, this prediction con- 
firms that dynamics will not be able to unambiguously 
reveal the presence of dark matter in most individual 
UCDs. Less luminous UCDs have less baryonic mass, 
and so will be more dynamically affected by the presence 
of dark matter if they reside in halos similar to those of 
more luminous UCDs. We also predict that among UCDs 
of similar luminosity, those with larger scale-sizes should 
have systematically higher dark matter fractions. This 
prediction makes sense, because larger half-light radii en- 
close a larger fraction of an object's dark matter halo, if 
UCDs of similar luminosity reside in similar dark matter 
halos. Current observations do not bear a clear signature 
of this predicted relationship (Mieske et al. 2008). How- 
ever, because of possible system-to-system variations and 
uncertainties in stellar M/L, it is impossible (to date) to 
draw robust conclusions about the dynamical evidence 
for dark matter or lack thereof. 

In making the quantitative predictions in Figure 2, we 
have assumed that UCDs contain the same amount of 
dark matter within their half-light radii as known dwarf 
galaxies with the same half-light radii. The dark matter 
halos inhabited by UCDs may instead have higher mass 
density than those inhabited by MW ultra-faint dwarfs, 
owing to gravitational contraction. Alternatively they 
may have lower mass density, owing to the far greater 
amount of feedback from star formation and death ex- 
perienced by UCDs with orders of magnitude more stars 
than ultra-faint dwarfs. Nevertheless, Figure 2 demon- 
strates a reasonable model in which dark matter is not 
dynamically detectable in most UCDs, but may be de- 
tectable in the least luminous UCDs. The relationship 
we predict between half-light radius and dynamical mass 
is dependent only on the assumption that similar lumi- 
nosity UCDs inhabit similar dark matter halos. 

4.2.2. UCD [Fe/H] 

Even if it is possible to assess crpe/H] in UCDs with 
Mv < —10, it would not easily a id in their classifica- 
tion (see §3.2). IBrodie et al.l (:2011') have recently argued 
that objects with lower stellar masses are also part of 
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Fig. 2. — The predicted {M/L)f^alf of UCDs with rf^alf = 30 or 75 pc, assuming they reside in dark matter halos like those inferred for 
Segue 1 and ComBer, respectively. Typical UCDs should not display dynamical evidence for dark matter, even if they do reside in the 
centers of dark matter halos. 

around the MW and other galaxies. For example, the 
spatial distribution of MW halo GCs is consistent with 
the predicted prese nt-day distribution of early forming 
dark matter peaks (jBrodie fc Straderll2006t iMoore et alj 
'20065. This similarity could be interpreted as evidence 
that GCs themselves reside in the center of present day 
dark matter halos and, if so, should be included in stud- 
ies that rely on dwarf galaxies as luminous tracers of the 
spatial and mass distribution of dark matter. 

No dynamical study of GCs has yielded a dynamical 
mass in excess of s tellar mass, even for l ower surface den- 
sity (Palomar 13, iBradford et al. '20 11D and tidally dis- 
rupting clusters (Palomar 5. Odenki rchen et al.l2003 ). In 
light of the dynamical arguments presented for UCDs, 
GCs would be unlikely to exhibit straightforward dy- 
namical evidence for dark matter even if they did reside 
in dark matter halos. The [Fe/H] analysis in §3.2 and 
shown in Figure 1 instead provides direct evidence that 
GCs do not satisfy the definition of a galaxy — the iron 
abundances of their stars is explicable with only stellar 
mass and Newtonian gravity. 

Additional indirect diagnostics also demonstrate that 
GCs would be classified as star clusters with our pro- 
posed definition. The pr esence of tid al streams around 
numerous MW GCs (e.g. ILeon et all [2000 ) provides up- 
per limits to their present-day masses; this is addi- 
tional evidence that their present-day dynamics are con- 
sistent with their observed stellar masses and Newto- 
nian gravity. Another diagnostic is the existence of GCs 
in low-mass dwarf galaxies, such as the Fornax dwarf 
spheroidal. If its GCs were embedded in dark matter ha- 
los, then their dynamical friction timescale for destruc- 
tion would be <1 Gyr, far shorter than their observed 



the UCD population around M87 in Virgo. NGC 2419, 
a MW GC, has a size (21 pc, Harris 1996) and absolute 
magnitude (My = -9.42. (Harris, 19 96) consistent with 
the lower luminosity UCDs around M87. At face value, 
NGC 2419's lack of an [Fe/H] spread (Table 1) suggests 
that star clusters may form with the sizes and luminosi- 
ties of at least some UCDs. However, NGC 2419's spread 
in Ca (^^0.2 dex) may be difficult to reconcile with the 
inferred depth of its potential well. Unlike spreads in 
lighter element s, a Ca spread mig ht require enrichment 
by supernovae ()Cohen et al.ll201(J ). 

It would be extremely interesting if future studies 
could measure (or set limits on) the [Fe/H] spread 
of a set of lower luminosity UCDs to see whether 
they all lack a spread in [Fe/H], as observed for 
typical star clusters. Another hint to a possible 
UCD-dwarf galaxy connection — or lack thereof — may be 
their average [Fe/H]. UCDs fall above the metallicity- 
luminosity relationship followed by dwarf galaxies (e.g., 
iChilingarian et al"]l201lL see also discussion in §3.3). If 
the UCDs are stripped remnants of nucleated dwarfs then 
they once would have been more luminous and may have 
fallen on observed metallicity-luminosity relationships. 

4.3. Globular Clusters 

A combination of dynamics, (7[fo/h]i and several in- 
direct diagnostics show that GCs, as a population, do 
not satisfy our definition of galaxy and do not presently 
inhabit dark matter halos. We briefly discuss this evi- 
dence here, because we should neither take for granted 
that canonical GCs do not satisfy our proposed defini- 
tion of galaxy, nor take for granted that they should 
be ignored in efforts to map dark matter substructure 
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ages (|Conrov fc Spergeill2011[ ). One final diagnostic may 
be the outer density p rofiles of GCs, as demonstrated 
by IConrov et all (|20T1 for the case of NGC 2419 and 
MGCl. 

Light element abundance spreads are common in GCs, 
and usually attributed to enrichment by asymptotic gi- 
ant b ranch stars or the winds of rotating mass ive stars 
(e.g.. iRenzinll [20081: Ne^tma fc D'AntoMll2009t) . These 
ejecta are less energetic than those of supernovae and can 
be retained by the gravity of stars alone. The ubiquity 
of these abundance variations, often identified through 
the anti-correlation of Na and O, has led to the sug- 
gestio n that such variations should define the class of 
GCs (jCarrettaet al.ll2010bl ). We do not advocate for 
this suggestion, since little is known about the abun- 
dance patterns of low-mass GCs, which may differ from 
those of more massive clusters, and star clusters with 
masses < W^Mq in t he Large Magellani c Cloud do not 
appear to self-enrich (jMilone et al.l 120091 ). Furthermore, 
more massive objects that might be confused with GCs, 
such as UCDs or dwarf nuclei, lack detailed abundance 
observations. 

Existing diagnostics do not preclude the hypothesis 
that some massive (My < -10) GCs may reside in 
dark matter halos. This possibility must be considered 
when comparing observations against cosmological mod- 
els. Extended star clusters (My ~ —7 to —8, rhaif ~ 20 - 
30 pc, Tanvir et al. 2012), such as those observed around 
M31 (JHuxor et al.1 12005'). also present a challenge to clas- 
sification. The M31 extended GCs would make particu- 
larly interesting targets for spectroscopic [Fe/H] studies, 
because their current stellar masses and escape velocities 
are too low to expect self-enrichment in iron. 

4.4. Tidal Dwarfs 

The term "tidal dwarf" (TD) refers to a gravitationally 
bound, galaxy-sized object (few kpc scale) fo rmed as a 
resul t of the tidal interaction of two galaxies ([Bournaudl 
l2010f) . These objects form from a combination of star 
formation in gaseous tidal tails and of the agglomera- 
tion of_existing_stars from the interacting parent galax- 
ies ([Kavirai et al.l [20121 . Candidates for such objects 
were originally o bserved in the Anten n ae and in compact 
galax y groups (jMirabel et al.l 119921 : iHunsberger et al.l 
|1996[) . Although many candidate TDs have been discov- 
ered since then, it remains difficult to determine whether 
TD candidates are truly self-bound (Due ct al. 20(13). 

Dynamical studies of TDs do not provide a defini- 
tive classification of these objects. Their kinematic 
properties are difficult to study, in part because TDs 
are typically observed while still embedded in ambi- 
ent tidal material from which they formed/are form- 
ing. Some studies find their dynamical masses to 
be consistent with their stellar and gas ( both neutral 
and niolecular) contents (iDuc et aLJ l2000t iBraine et al.l 
[200ll:lBmiFnaud et al.l 12004 iDuc et al.ll2007[ ). while oth- 
ers find dynamical masses 2 — 3 times higher than ex- 
pecte d from observed stars and gas (jBournaud et al.l 
I2007f ). In all cases, the uncertainties are substantial. 
Even in a cold dark matter interpretation of galaxies, 
TDs are not expected con tain (much) dark matter (e.g., 
iBarnes fc HernquTs^ll992D . Unlike gas, the dark mat- 
ter in TD progenitor material cannot dissipate energy 
and has a velocity dispersion exceeding the escape ve- 



locity of the forming TD ()Bournaudll2010f) . unless some 
dark matter i s pres e nt in a cold, rota ting, galaxy disk 
(jPurceU et al.l l2009t iRead et al.ll2009D . Identifying a 
sample of relatively older (>1 Gyr) TDs and conduct- 
ing uniform dynamical studies will help reveal whether: 
(i) TDs are simply composed of gas and stars orbiting in 
a Newtownian potential, (ii) galaxy disks do contain a 
dark matter component which can be accreted by form- 
ing TDs, or (iii) TDs demonstrate a dynamical regime 
governed by non-Newtonian gravity. If (ii) or (iii) is ver- 
ified, then TDs would be classified by galaxies by our 
definition. 

The possible contribution of ancient TDs formed at 
high redshift to today's dwarf galaxy population, in par- 
ticular around the Milky Way, is controversial. Observa- 
tions of the universe at low or intermediate redshift imply 
that that TDs could not contribute more tha n ^10% of 
the dwarf galaxies in the local universe (e.g.. lWen et al.l 
120121: iKavirai et al.|[20T2l ). TDs forming in the local uni- 
verse also do not exhibit the relations hip between stel- 
lar mass and metallicity (iWeilbacher et al. 2003) that is 
observed in the MW dwarfs (|Kirbv et al.l 120111 ). More- 
over, kinematic studies of nearby TDs do not imply the 
high dynamical M/L observed for MW dwarf satellites. 
Others propose that the MW's dwarf galaxies may be 
dominated by tidal dwarfs formed at very high redshift 
when merger rates were far higher, and that the high 
M/L inferred for MW dwarfs are actuall y a misinterpre- 
tation of the observed kine matics (e.g. Kuh nfc Milled 
[19891: [Meti fc Kroupall2"007t iKroupa et al. 20i0l ). It re- 
mains to be seen whether models of ancient TDs evolving 
into z = dwarfs could fall on the same metallicity- 
luminosity relation followed by both MW dwarfs and 
spheroidal galaxies over a wide range of masses. 

5. DISCUSSION AND CONCLUSIONS 

To compound the ambiguities inherent to classify- 
ing objects such as extreme ultra-faint MW dwarfs and 
UCDs, observers have neither agreed upon a definition 
of galaxy nor reached consensus on how to interpret ob- 
servations in hand. To facilitate comparisons between 
dwarf galaxy predictions and the increasingly complex 
sets of observations of candidate dwarf galaxies, the field 
needs an agreed-upon definition for galaxy. We have ac- 
cordingly proposed a physically motivated definition that 
does not insist on a cold dark matter interpretation of 
data: A galaxy is a gravitationally bound collection of 
stars whose properties cannot be explained by a combina- 
tion of baryons and Newton's laws of gravity. 

We have explored possible diagnostics of this galaxy 
definition (primarily in the context of a cold dark mat- 
ter dominated universe) , primarily kinematic studies and 
[Fe/H] spread. Although kinematic studies generally pro- 
vide the most direct way to infer a galaxy definition, it 
can be difficult to measure the dynamical mass of low lu- 
minosity and/or low velocity dispersion (< few km s"""^) 
systems. Even once robustly measured, interpreting rel- 
atively modest dynamical M/L (< 10) may face several 
stumbling blocks: those that could have generated over- 
estimates of dynamical mass (e.g., binary stars, contam- 
inants in spectroscopic sample) and those could gener- 
ate underestimates of stellar mass (e.g., sparse sampling 
of the stellar luminosity function, an overabundance of 
stellar remnants). While these effects do not appear to 
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be a major problem for objects currently classified as 
galaxies, including the Milky Way's dwarfs, they should 
be carefully considered as discoveries at the extremes of 
the cosmic zoo continue. Systems such as UCDs and 
massive GCs also may not bear a kinematic signature of 
dark matter or non-Newtonian gravity, even if present, 
because their baryons are so densely packed. 

C[Fo/H] provides an complimentary means to diagnose 
a galaxy definition for systems less luminous than My = 
— 10. Using public spectroscopic [Fe/H] measurements, 
we recalculated the average systemic [Fe/H] and associ- 
ated dispersions for 24 Milky Way GCs and 16 Milky 
Way dwarf galaxies. All dwarf galaxies show spectro- 
scopic [Fe/H] spreads of ~0.3 dex or more. No GC less 
luminous than My — —10 shows a notable (> 0.1 dex) 
[Fe/H] spread. The a\F c/R] diagnost ic has already been 
applied to the Seg ue 1 (iSimon et al.ll2011h and Willman 
1 dwarf galaxies (|Willman et al.l 1201 Ih . One possible 
caveat with the (T[Fe/H] diagnostic is the possibility that 
the mergers of multiple star clusters could yield an iron 
abundance spread. This merging star cluster hypothesis, 
which would produce a multimodal [Fe/H] distribution, 
should be carefully considered when classifying objects 
by [Fe/H] dispersion alone. 

The Fundamental Plane and its variants do not 
presently provide an alternative means to diagnose a 
galaxy definition for low luminosity systems. However, 
these scaling relations do provide a useful benchmark 
against which to compare ambiguous objects. For ex- 
ample, an outlier from known scaling relations may sig- 
nal a problem with its calculated velocity dispersion or 
estimated stellar mass (such as the issues discussed in 
§3.1.1.) Well behaved scaling relations can also help 
rule out pathological explanations for sets of objects, 
especially when metallicity is included. For example, 
the metallicity-luminosity relation followed by the Milky 
Way's lowest luminosity dwarfs helps rules out alterna- 
tive hypotheses for their existence as a population, such 
as tidal tails at apocenter and clumps in streams. 

There are some classes of objects not discussed in this 
paper, but which would be worth new consideration in 
the context of our proposed galaxy definition. For exam- 
ple, dwarf ellipticals (dEs) do not typically show strong 
kinematic evidenc e for non-baryonic ma ss in their cen - 
tral regions (e.g. iWolf et al. 2010; For bes et al.l 1201 ID . 
However, recent kinematic studies of stars and globu- 
lar clusters in th eir outer regions (jBeaslev et al.l 120091 : 
iGeha et al]l2010[ ) have consistently suggested that M/L 
increases with radius and that stars alone cannot account 
for the observations. Therefore, the current data favors 
the classification of dEs as galaxies by our definition. If 
future observations of dEs do not support this emerging 
consensus, then this classification should be revisited. 

After examining massive globulars, UCDs and tidal 
dwarfs in detail, we find that they can not yet con- 
clusively be classified given existing diagnostics of our 
galaxy definition. Their ultimate classification must be 
guided by future observational data. If UCDs and tidal 
dwarfs are inconsistent with a galaxy definition, this does 
not mean that they should automatically be classified as 
star clusters. Both of these classes of objects are interest- 
ing and stand on their own as worthy to investigate, given 
their unique properties and possibly formation channels. 



For tidal dwarfs, in particular, we advocate that these 
objects are not lumped in with galaxies or clusters but 
that they remain their own distinct class of objects. 

We have suggested several measurements, some of 
which are possible now, that could facilitate the classifi- 
cation of these, and other, extreme objects: 

• Observational constraints on the [Fe/H] spread in 
extended GCs and any UCD fainter than My = —10 
(§3.2.2, §4.3); 

• Adaptive optics spectroscopy to measure the [Fe/H] of 
individual stars in massive M31 GCs (§3.2.3, §4.3); 

• Dynamical studies of the lowest luminosity (My > —9), 
largest scale sized UCDs (vhaif ^ 30 pc) (§4.2); 

• Measuring dynamical masses of a larger sample of 
UCDs to look for a positive correlation between half-light 
radius and half-light M/L, at set UCD luminosity (§4.2); 

• Dynamical and chemical studies of extreme MW 
satellites Bootes II and Segue 2 (§4.1). 

Basic cold dark matter plus galaxy formation models 
predict a dichotomy between systems that form in the 
centers of dark matter halos and systems that form in the 
monolithic collapse of gas clouds that are not the primary 
baryonic components of dark matter halos. Our best 
present interpretation of the observations in this con- 
text reveals that systems forming at the center of a dark 
matter halo bear an observable imprint of this formation 
channel, such as the kinematic or chemical diagnostics 
discussed here. This observable imprint would translate 
to a galaxy classification by our proposed definition. The 
fact that systems classified as galaxies may be equivalent 
to the set of astrophysical systems that formed in dark 
matter halos can be used as strong guidance to theo- 
rists when selecting against which astrophysical systems 
to compare their predictions. However, even in a dark 
matter context, we cannot take for granted that systems 
classified as galaxies by our definition are inclusive of all 
systems that formed in dark matter halos and exclusive 
of systems that formed otherwise. 

As our understanding of the universe grows, it may 
be possible for systems that formed inside of dark mat- 
ter halos to fail the galaxy diagnostics discussed here. 
For example, it would not be unreasonable to conceive 
that a very low-luminosity fossil galaxy could form all 
of its stars over a sufficiently short timescale that no 
opportunity for self-enrichment by supernovae could oc- 
cur, leading to a minimal spread in [Fe/H]. If this is the 
case, then alternative diagnostics need to be identified 
for such "first" galaxies. No objects meeting our defini- 
tion of a galaxy via kinematics, but without a spread in 
[Fe/H] , have yet been discovered, but it is plausible they 
exist. Cosmolog ical globular cluster s, if they exist, may 
be such objects (jGriffen et al.lf2010l ). Conversely, it may 
be possible for stellar systems that formed inside of a 
dark matter halo to lose most or all of their dark matter. 
In this scenario, a(n almost) stellar-only system may ex- 
ist with the chemical imprint of formation within a dark 
matter halo. Although simulations of dwarf galaxies in 
both cuspy and cored halos show that this is unlikely (see 
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§3.1), simulations of globular clusters within dark mat- 
ter halos have shown that it may be possible to remove 
most of the dark matter in sy stems close to disruption 
(e.g. IMashchenko fc SillsllMiHh . 

Our proposed galaxy definition is itself independent of 
our observational knowledge and currently favored the- 
ories for structure formation; it can thus remain un- 
changed even as our understanding of the complex uni- 
verse evolves. However, the particular diagnostics of 
this definition as investigated in this polemic may indeed 
need to be revisited as our knowledge of extreme objects 
grows — both observationally and theoretically. For ex- 
ample, the possible use of spreads in elements other than 
iron (such as calcium) to diagnose a galaxy classification 
is something that should continue to be scrutinized as 
our knowledge of such abundance patterns grow. 
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TABLE 1 

[Fe/H] properties of MW globular clusters and dwarfs 



Name 



[Fe/H] 
dcx 



±34% CL 
dex 



'^[Fo/H] 
dex 



+34% CL 
dex 



-34% CL 
dcx 



Mv 



Ref 



type 



ijCen 


-1.647 


0.009 


0.271=' 


0.007 


0.007 


-10.3 


855 


JIO 


GC 


M54 


-1.559 


0.021 


0.186 


0.016 


0.014 


-10.0 


76 


Carlo 


GC 


NGC 6441 


-0.334 


0.018 


0.079 


0.016 


0.013 


-9.6 


25 


G07 


GC 


NGC 104 


-0.743 


0.003 


0.024 


0.003 


0.002 


-9.4 


147 


Car09b 


GC 


NGC 2419 


-2.095 


0.019 


0.032 


0.013 


0.009 


-9.4 


38 


CohlO 


GC 


NGC 2808 


-1.105 


0.006 


0.062 


0.005 


0.004 


-9.4 


123 


Car06 


GC 


NGC 6388 


-0.404 


0.014 


0.071 


0.012 


0.010 


-9.4 


36 


Car09b 


GC 


NGC 7078 


-2.341 


0.007 


0.055 


0.006 


0.005 


-9.2 


84 


Car09b 


GC 


NGC 5904 


-1.346 


0.002 


0.014 


0.002 


0.002 


-8.8 


136 


Car09b 


GC 


M22 


-1.764 


0.016 


0.099*= 


0.013 


0.011 


-8.5 


37 


Mil 


GC 


NGC 1851 


-1.157 


0.005 


0.046 


0.004 


0.003 


-8.3 


124 


Carll 


GC 


NGC 1904 


-1.545 


0.005 


0.028 


0.005 


0.004 


-7.9 


58 


Car09b 


GC 


NGC 6752 


-1.564 


0.004 


0.034 


0.003 


0.003 


—7.7 


137 


Car07b 


GC 


NGC 6809 


-1.970 


0.004 


0.035 


0.003 


0.003 


-7.6 


156 


Car09b 


GC 


NGC 3201 


-1.495 


0.004 


0.042 


0.004 


0.003 


-7.5 


149 


Car09b 


GC 


NGC 6254 


-1.557 


0.005 


0.048 


0.004 


0.003 


-7.5 


147 


Car09b 


GC 


NGC 7099 


-2.358 


0.006 


0.037 


0.006 


0.005 


-7.5 


65 


Car09b 


GC 


NGC 4590 


-2.230 


0.007 


0.057 


0.006 


0.005 


-7.4 


122 


Car09b 


GC 


NGC 6218 


-1.313 


0.004 


0.027 


0.004 


0.003 


-7.3 


79 


Car07a 


GC 


NGC 6121 


-1.200 


0.003 


0.018 


0.003 


0.002 


-7.2 


103 


Car09b 


GC 


NGC 6171 


-1.066 


0.008 


0.037 


0.007 


0.006 


-7.1 


33 


Car09b 


GC 


NGC 288 


-1.219 


0.004 


0.034 


0.004 


0.003 


-6.8 


110 


Car09b 


GC 


NGC 6397 


-1.994 


0.004 


0.027 


0.003 


0.003 


-6.6 


144 


Car09b 


GC 


NGC 6838 


-0.806 


0.006 


0.027 


0.005 


0.005 


-5.6 


39 


Car09b 


GC 


For 


-1.025 


0.012 


0.292 


0.010 


0.010 


-13.3 


675 


KIO 


dwarf 


Leo I 


-1.450 


0.011 


0.276 


0.009 


0.008 


-11.9 


827 


KIO 


dwarf 


Scl 


-1.726 


0.024 


0.452 


0.019 


0.017 


-11.2 


376 


KIO 


dwarf 


Leo II 


-1.670 


0.024 


0.347 


0.020 


0.018 


-10.0 


258 


KIO 


dwarf 


Sex 


-1.966 


0.039 


0.339 


0.033 


0.030 


-9.6 


141 


KIO 


dwarf 


Dra 


-1.946 


0.024 


0.354 


0.020 


0.019 


-8.8 


298 


KIO 


dwarf 


CVnl 


-1.962 


0.038 


0.441 


0.032 


0.029 


-8.6 


174 


KIO 


dwarf 


UMi 


-2.112 


0.027 


0.319 


0.025 


0.023 


-9.2 


212 


KIO 


dwarf 


Here 


-2.518 


0.140 


0.583 


0.124 


0.095 


-6.2 


21 


K08 


dwarf 


UMal 


-2.334 


0.128 


0.638 


0.106 


0.086 


-5.5 


31 


K08 


dwarf 


Leo IV 


-2.363 


0.230 


0.695 


0.210 


0.149 


-5.5 


12 


K08 


dwarf 


CvnII 


-2.444 


0.178 


0.621 


0.164 


0.120 


-4.6 


15 


K08 


dwarf 


UMall 


-2.357 


0.204 


0.563 


0.204 


0.136 


-4.0 


9 


K08 


dwarf 


ComBer 


-2.640 


0.100 


0.382 


0.088 


0.070 


-3.8 


23 


K08 


dwarf 


Will 


-2.110 


0.367 


0.557 


0.577 


0.231 


-2.7 


3 


Wll 


dwarf 


Segl 


-2.735 


0.389^^ 


0.752 


0.417 


0.227 


-1.5 


7 


NIO, Sll 


dwarf 



Note. — The reference column gives the source of individual [Fc/H] measurements used to estimate the 
dispersion in each object. For Segue 1, on l y the one star (Seg 1-7) is taken from Norris ct al. (201C|). Values 
of My for the dwarfs are from [Sand et all J2Q11D and references therein. Values of Mv for the GCs are from 

I Harris! l|T996;, 2010 edition). The posterior distribution of [Fe/H] sufficiently symmetric that we only quote a 
single value for ±34% CL, taking the average of the + and — values in t he small number of cases with a few 
thousandth of a dex d ifference between the two . . Refer ence key: JIO — [Johnson & Pilachowski (201 y), Carll — 

ICarretta et al.| ) |201lD. CarlO ^[Carretta et all | |2010aD, CohlO ^ICohen e t al. (2010), Car09b ^ Carretta et al.| 
I2p09b). Car06 ^ Carrctta et al. (2006), Mil ^ Marino et al. (2011), Car07a ^ Ca rretta"et~al| J2007all G07 ^ 

[Gratton et al. (2007), Car07b ^ Carret ta ct al.. (^20 07b), KIO ^ Kirby_ct al. ( 2 010,) , K08 ^ IKirby et al.| j200gD . 
Wll ^ iWillman et al.i (|2011|) . NIO ^ iNorris et al.h201Q) . Sll ^ iSimon et al.| l l201ll) 
^ This value is a lower limit (see §3.2.1). 

This value is an upper limit (see §3.2.1). 
^ Unlike the other objects, the metallicity of Segue 1 has asymmetric uncertainties: — 2.735_Q'4Qg 



14 



Willman & Strader 



REFERENCES 



Aihara, H. e. a. 2011, ApJS, 193, 29 

Bailin, J., & Harris, W. E. 2009, ApJ, 695, 1082 

Barnes, J. E., & Hernquist, L. 1992, Nature, 360, 715 

Baumgardt, H., Kroupa, P., & Parmentier, G. 2008, MNRAS, 

384, 1231 
Baumgardt, H., & Makino, J. 2003, MNRAS, 340, 227 
Beasley, M. A., Cenarro, A. J., Strader, J., & Brodie, J. P. 2009, 

AJ, 137, 5146 
Bekki, K., & Freeman, K. C. 2003, MNRAS, 346, Lll 
Belokurov, V. et al. 2009, MNRAS, 397, 1748 
Bournaud, F. 2010, Advances in Astronomy, 2010 
Bournaud, F., Due, P.-A., Amram, P., Combes, F., & Gach, J.-L. 

2004, A&A, 425, 813 
Bournaud, F. et al. 2007, Science, 316, 1166 
Boylan-Kolchin, M., Bullock, J. S., & Kaplinghat, M. 2011, 

MNRAS, L267+ 
— . 2012, MNRAS, 422, 1203 
Bradford, J. D. et al. 2011, ApJ, 743, 167 
Braine, J., Due, P.-A., Lisenfeld, U., Charmandaris, V., Vallejo, 

O., Leon, S., & Brinks, E. 2001, A&A, 378, 51 
Brodie, J. P., Romanowsky, A. J., Strader, J., & Forbes, D. A. 

2011, AJ, 142, 199 
Brodie, J. P., & Strader, J. 2006, ARA&A, 44, 193 
Busha, R., Lake, G., & Reed, D. 2012, in preparation 
Carretta, E., Bragaglia, A., Gratton, R., D'Orazi, V., & 

Lucatello, S. 2009a, A&A, 508, 695 
Carretta, E. et al. 2007a, A&A, 464, 939 
Carretta, E., Bragaglia, A., Gratton, R. G., Leone, F., 

Recio-Blanco, A., & Lucatello, S. 2006, A&A, 450, 523 
Carretta, E. et al. 2010a, A&A, 520, A95 
— . 2009b, A&A, 505, 117 
Carretta, E., Bragaglia, A., Gratton, R. G., Lucatello, S., & 

Momany, Y. 2007b, A&A, 464, 927 
Carretta, E., Bragaglia, A., Gratton, R. G., Recio-Blanco, A., 

Lucatello, S., D'Orazi, V., & Cassisi, S. 2010b, A&A, 516, A55 
Carretta, E., Lucatello, S., Gratton, R. G., Bragaglia, A., & 

D'Orazi, V. 2011, A&A, 533, A69 
Chiboucas, K. et al. 2011, ApJ, 737, 86 
Chilingarian, L V., Mieske, S., Hilker, M., & Infante, L. 2011, 

MNRAS, 412, 1627 
Cohen, J. G., Kirby, E. N., Simon, J. D., & Geha, M. 2010, ApJ, 

725, 288 
Cohen, J. G., & Melendez, J. 2005, AJ, 129, 303 
Conroy, C, Loeb, A., & Spergel, D. N. 2011, ApJ, 741, 72 
Conroy, C, & Spergel, D. N. 2011, ApJ, 726, 36 
Crommelin, A. C. D. 1918, JRASC, 12, 33 
Da Rocha, C, Mieske, S., Georgiev, I. Y., Hilker, M., Ziegler, 

B. L., & Mendes de Oliveira, C. 2011, A&A, 525, A86 
Dabringhausen, J., Hilker, M., & Kroupa, P. 2008, MNRAS, 386, 

864 
Dabringhausen, J., Kroupa, P., & Baumgardt, H. 2009, MNRAS, 

394, 1529 
Dabringhausen, J., Kroupa, P., Pflamm-Altenburg, J., & Mieske, 

S. 2012, ApJ, 747, 72 
D'Antona, F., Bellazzini, M., Caloi, V., Pecci, F. F., Galleti, S., & 

Rood, R. T. 2005, ApJ, 631, 868 
de Sitter, W. 1917, MNRAS, 78, 3 
Dehnen, W., Odenkirchen, M., Grebel, E. K., & Rix, H. 2004, AJ, 

127, 2753 
Dekel, A., Stoehr, F., Mamon, G. A., Cox, T. J., Novak, G. S., & 

Primack, J. R. 2005, Nature, 437, 707 
Dopita, M. A., & Smith, G. H. 1986, ApJ, 304, 283 
Drinkwater, M. J., Gregg, M. D., Hilker, M., Bekki, K., Couch, 

W. J., Ferguson, H. C, Jones, J. B., & Phillipps, S. 2003, 

Nature, 423, 519 
Due, P.-A. 2012, Birth, Life and Survival of Tidal Dwarf 

Galaxies, in Dwarf Galaxies: Keys to Galaxy Formation and 

Evolution, Astrophysics and Space Science Proceedings, ISBN 

978-3-642-22017-3. Springer- Verlag Berlin Heidelberg, 2012, 

p. 305 (eds. Papaderos, P. and Recchi, S. and Hensler, G.) 
Due, P.-A., & Bournaud, F. 2008, ApJ, 673, 787 
Due, P.-A., Braine, J., Lisenfeld, U., Brinks, E., & Boquien, M. 

2007, A&A, 475, 187 
Due, P.-A., Brinks, E., Springel, V., Pichardo, B., Weilbacher, P., 

& Mirabel, I. F. 2000, AJ, 120, 1238 



Evstigneeva, E. A., Gregg, M. D., Drinkwater, M. J., & Hilker, 

M. 2007, AJ, 133, 1722 
Fadely, R., Willman, B., Geha, M., Walsh, S., Murioz, R. R., 

Jerjen, H., Vargas, L. C, & Da Costa, G. S. 2011, AJ, 142, 88 
Fellhauer, M., & Kroupa, P. 2005, MNRAS, 359, 223 
Ferraro, F. R. et al. 2009, Nature, 462, 483 
Firth, P., Evstigneeva, E. A., & Drinkwater, M. J. 2009, MNRAS, 

394, 1801 
Forbes, D. A., & Kroupa, P. 2011, PubUcations of the 

Astronomical Society of Australia, 28, 77 
Forbes, D. A., Lasky, P., Graham, A. W., & Spitler, L. 2008, 

MNRAS, 389, 1924 
Forbes, D. A., Spitler, L. R., Graham, A. W., Foster, C, Hau, 

G. K. T., & Benson, A. 2011, MNRAS, 413, 2665 
Forbes, D. A., Spitler, L. R., Harris, W. E., Bailin, J., Strader, J., 

Brodie, J. P., & Larsen, S. S. 2010, MNRAS, 403, 429 
Frank, M. J., Hilker, M., Mieske, S., Baumgardt, H., Grebel, 

E. K., & Infante, L. 2011, MNRAS, 414, L70 
Fuentes-Carrera, I., Jablonka, P., Sarajedini, A., Bridges, T., 

Djorgovski, G., & Meylan, G. 2008, A&A, 483, 769 
Geha, M., van der Marel, R. P., Guhathakurta, P., Gilbert, 

K. M., Kalirai, J., & Kirby, E. N. 2010, ApJ, 711, 361 
Geha, M., Willman, B., Simon, J. D., Strigari, L. E., Kirby, 

E. N., Law, D. R., & Strader, J. 2009, ApJ, 692, 1464 
Giersz, M. 2001, MNRAS, 324, 218 
Gilmore, G., Wilkinson, M. I., Wyse, R. F. G., Kleyna, J. T., 

Koch, A., Evans, N. W., & Grebel, E. K. 2007, ApJ, 663, 948 
Gnedin, O. Y., Zhao, H., Pringle, J. E., Fall, S. M., Livio, M., & 

Meylan, G. 2002, ApJ, 568, L23 
Governato, F. et al. 2010, Nature, 463, 203 
Graham, A. W. 2011, ArXiv e-prints 1108.0997, condensed 

version of a review to appear in "Planets, Stars and Stellar 

Systems", Springer pub. 2012 
Gratton, R., Sneden, C, & Carretta, E. 2004, ARA&A, 42, 385 
Gratton, R. G. et al. 2007, A&A, 464, 953 
Griffen, B. F., Drinkwater, M. J., Thomas, P. A., Helly, J. C, & 

Pimbblet, K. A. 2010, MNRAS, 405, 375 
Ha§egan, M. et al. 2005, ApJ, 627, 203 
Harris, W. E. 1996, AJ, 112, 1487 
Harris, W. E., Whitmore, B. C, Karakla, D., Okoh, W., Baum, 

W. A., Hanes, D. A., & Kaveiaars, J. J. 2006, ApJ, 636, 90 
Hernandez, X. 2012, MNRAS, 420, 1183 
Hilker, M., Infante, L., Vieira, G., Kissler-Patig, M., & Richtler, 

T. 1999, A&AS, 134, 75 
Hubble, E. P. 1926, ApJ, 64, 321 
Hunsberger, S. D., Charlton, J. C, & Zaritsky, D. 1996, ApJ, 

462, 50 
Huxor, A. P., Tanvir, N. R., Irwin, M. J., Ibata, R., CoUett, J. L., 

Ferguson, A. M. N., Bridges, T., & Lewis, G. F. 2005, MNRAS, 

360, 1007 
Ivezie, Z., Tyson, J. A., AUsman, R., Andrew, J., A ngel, R., & for 

the LSST Collaboration. 2008, apCiv:0805.2366l also available 

at http://www.lsst.org/overviewl 
Johnson, C. 1., & Pilaehowski, C. A. 2010, ApJ, 722, 1373 
Kaiser, N. ct al. 2002, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 4836, 

Society of Photo-Optical Instrumentation Engineers (SPIE) 

Conference Series, ed. J. A. Tyson & S. Wolff, 154-164 
Kaviraj, S., Darg, D., Lintott, C, Schawinski, K., & Silk, J. 2012, 

MNRAS, 419, 70 
Keller, S. C. et al. 2007, Publications of the Astronomical Society 

of AustraUa, 24, 1 
Kirby, E. N. et al. 2010, ApJS, 191, 352 
Kirby, E. N., Lanfranchi, G. A., Simon, J. D., Cohen, J. G., & 

Guhathakurta, P. 2011, ApJ, 727, 78 
Kirby, E. N., Simon, J. D., Geha, M., Guhathakurta, P., & 

Frebel, A. 2008, ApJ, 685, L43 
Klimentowski, J., Lokas, E. L., Kazantzidis, S., Prada, F., Mayer, 

L., & Mamon, G. A. 2007, MNRAS, 378, 353 
Koposov, S. E. et al. 2011, ApJ, 736, 146 
Kroupa, P. 2008, in lAU Symposium, Vol. 246, lAU Symposium, 

ed. E. Vesperini, M. Giersz, & A. Sills, 13-22 
Kroupa, P. et al. 2010, A&A, 523, A32 
Kuhn, J. R., & Miller, R. H. 1989, ApJ, 341, L41 



■'Galaxy," defined 



15 



Lee, Y., Joo, J., Sohn, Y., Rey, S., Lee, H., & Walker, A. R. 1999, 

Nature, 402, 55 
Leon, S., Meylan, G., & Combes, F. 2000, A&A, 359, 907 
Madrid, J. P. et al. 2010, ApJ, 722, 1707 
Majewski, S. R., Nidever, D. L., Smith, V. V., Damke, G. J., 

Kunkel, W. E., Patterson, R. J., Bizyaev, D., & Garcia Perez, 

A. E. 2012, ApJ, 747, L37 
Maraston, C., Bastian, N., Saglia, R. P., Kissler-Patig, M., 

Schweizer, F., &; Goudfrooij, P. 2004, A&A, 416, 467 
Marino, A. F. et al. 2011, A&A, 532, A8 
Martin, N. F., de Jong, J. T. A., & Rix, H.-W. 2008, ApJ, 684, 

1075 
Martin, N. F., Ibata, R. A., Chapman, S. C, Irwin, M., & Lewis, 

G. F. 2007, MNRAS, 380, 281 
Martinez, G. D., Minor, Q. E., Bullock, J., Kaplinghat, M., 

Simon, J. D., & Geha, M. 2011, ApJ, 738, 55 
Mashchenko, S., & Sills, A. 2005, ApJ, 619, 258 
McConnachie, A. W., & Cote, P. 2010, ApJ, 722, L209 
Messier, C. 1781, Connoissance des Temps for 1784 
Metz, M., & Kroupa, P. 2007, MNRAS, 376, 387 
Meylan, G., Sarajedini, A., Jablonka, P., Djorgovski, S. G., 

Bridges, T., & Rich, R. M. 2001, AJ, 122, 830 
Mieske, S. et al. 2008, A&A, 487, 921 
— . 2006, ApJ, 653, 193 
— . 2010, ApJ, 710, 1672 
Milgrom, M. 1983, ApJ, 270, 371 
Milone, A. P., Bedin, L. R., Piotto, G., & Anderson, J. 2009, 

A&A, 497, 755 
Minchin, R. et al. 2005, ApJ, 622, L21 

Mirabel, L F., Dottori, H., & Lutz, D. 1992, A&A, 256, L19 
Misgeld, I., & Hilker, M. 2011, MNRAS, 414, 3699 
Misgeld, I., Mieske, S., Hilker, M., Richtler, T., Georgiev, 1. Y., & 

Schuberth, Y. 2011, A&A, 531, A4 
Moore, B., Diemand, J., Madau, P., Zemp, M., & Stadel, J. 2006, 

MNRAS, 368, 563 
Muiioz, R. R., Geha, M., & Willman, B. 2010, AJ, 140, 138 
Norris, J. E., Wyse, R. F. G., Gilmore, G., Yong, D., Frebel, A., 

Wilkinson, M. L, Belokurov, V., & Zucker, D. B. 2010, ApJ, 

723, 1632 
Norris, M. A., & Kannappan, S. J. 2011, MNRAS, 414, 739 
Odenkirchen, M., Grebel, E. K., Dehnen, W., Rix, H.-W., & 

Cudworth, K. M. 2002, AJ, 124, 1497 
Origlia, L. et al. 2011, ApJ, 726, L20 
Penarrubia, J., Benson, A. J., Walker, M. G., Gilmore, G., 

McConnachie, A. W., & Mayer, L. 2010, MNRAS, 406, 1290 
Penarrubia, J., Navarro, J. F., & McConnachie, A. W. 2008, ApJ, 

673, 226 
Pontzen, A., & Governato, F. 2012, MNRAS, 421, 3464 
Purcell, C. W., Bullock, J. S., & Kaplinghat, M. 2009, ApJ, 703, 

2275 
Read, J. I., Mayer, L., Brooks, A. M., Governato, F., & Lake, G. 

2009, MNRAS, 397, 44 
Renzini, A. 2008, MNRAS, 391, 354 
Romanowsky, A. J., Douglas, N. G., Arnaboldi, M., Kuijken, K., 

Merrifield, M. R., Napolitano, N. R., Capaccioli, M., & 

Freeman, K. C. 2003, Science, 301, 1696 



Sand, D. J., Strader, J., Willman, B., Zaritsky, D., McLeod, B., 

CaldweU, N., Seth, A., & Olszewski, E. 2011, ArXiv e-prints 

1111.6608 
Sarajedini, A., & Layden, A. C. 1995, AJ, 109, 1086 
Saviane, I., Da Costa, G. S., Held, E. V., Sommariva, V., 

GuUieuszik, M., Barbuy, B., & Ortolani, S. 2012, ArXiv e-prints 
Shapley, H. 1919, PASP, 31, 261 
Simon, J. D., & Geha, M. 2007, ApJ, 670, 313 
Simon, J. D. et al. 2011, ApJ, 733, 46 
Skillman, E. D., Kennicutt, R. C, & Hodge, P. W. 1989, ApJ, 

347, 875 
Slipher, V. M. 1917, Proceedings of the American Philosophical 

Society, 56, 403 
Sotiriou, T. P., & Faraoni, V. 2010, Reviews of Modern Physics, 

82, 451 
Spitler, L. R., Larsen, S. S., Strader, J., Brodie, J. P., Forbes, 

D. A., & Beasley, M. A. 2006, AJ, 132, 1593 
Strader, J., Brodie, J. P., Spitler, L., & Beasley, M. A. 2006, AJ, 

132, 2333 
Strader, J., Caldwell, N., & Seth, A. C. 2011, AJ, 142, 8 
Strader, J., & Smith, G. H. 2008, AJ, 136, 1828 
Tanvir, N. R. et al. 2012, MNRAS, 422, 162 
The Dark Energy Survey Collaboration. 2005, White Paper 

submitted to the Dark Energy Task Force, 

larXiv:astro-ph/0510346 
Tollerud, E. J., Bullock, J. S., Graves, G. J., & Wolf, J. 2011, 

ApJ, 726, 108 
Tremonti, C. A. et al. 2004, ApJ, 613, 898 
van den Bergh, S. 2008, MNRAS, 385, L20 
Ventura, P., & D'Antona, F. 2009, A&A, 499, 835 
Vesperini, E., & Heggie, D. C. 1997, MNRAS, 289, 898 
Walker, M. G., Mateo, M., Olszewski, E. W., Penarrubia, J., 

Wyn Evans, N., & Gilmore, G. 2009a, ApJ, 704, 1274 
Walker, M. G., Mateo, M., Olszewski, E. W., Sen, B., & 

Woodroofe, M. 2009b, AJ, 137, 3109 
Walsh, S. M., Willman, B., Sand, D., Harris, J., Seth, A., 

Zaritsky, D., & Jerjen, H. 2008, ApJ, 688, 245 
Weilbacher, P. M., Due, P.-A., & Fritze-v. Alvensleben, U. 2003, 

A&A, 397, 545 
Wen, Z.-Z., Zheng, X.-Z., Zhao, Y.-H., & Gao, Y. 2012, Ap&SS, 

337, 729 
Willman, B., Geha, M., Strader, J., Strigari, L. E., Simon, J. D., 

Kirby, E., Ho, N., & Warres, A. 2011, AJ, 142, 128 
Wolf, J., Martinez, G. D., Bullock, J. S., Kaplinghat, M., Geha, 

M., Muiioz, R. R., Simon, J. D., & Avedo, F. F. 2010, MNRAS, 

406, 1220 
Woo, J., Courteau, S., & Dekel, A. 2008, MNRAS, 390, 1453 
Zaritsky, D., Zabludoff, A. I., & Gonzalez, A. H. 2011, ApJ, 727, 

116 
Zolotov, A., Hogg, D. W., & Willman, B. 2011, ApJ, 727, L14 



